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A selection of topics on open heavy quark production at HERA are reviewed here. Measure- 
ments of charm fragmentation parameters will be presented together with developments in 
the study of D* di-jet angular distributions. Charm production in deep inelastic scattering 
(DIS) is also discussed. Finally we deal with recent measurements of b cross sections using 
impact parameters in both DIS and photoproduction regimes. 



1 Introduction 

Colliding ep at a center of mass energy of yfs = 296 — 318 GeV HERA provides an interesting 
environment for testing QCD predictions on heavy quark production. The virtual photon emit- 
ted by the incoming lepton provides a clean probe which, interacting with quarks and gluons 
in the proton, can initiate hard processes. The scale of the QCD interaction spans over a wide 
range of values which are under direct experimental control. This report will concentrate on 
a selection of some recent measurements performed by the HI and ZEUS collaborations with 
a focus in particular on open c and b quark production in both deep inelastic scattering (DIS: 
Q 2 > 1 GeV 2 ) and photoproduction (Q 2 ~ GeV 2 ) regimes. 



2 Charm production 

2.1 Fragmentation tests 

The luminosity accumulated in the first phase of HERA running (1992-00, ~ 130 pb" 1 ) allows 
the study of particular decay channels which provide measurements of some phenomenological 
parameters used to describe charm fragmentation. ZEUS recently presented measurements of 
the branching fraction f(c — > D^) Q 0, the strangeness suppression factor i 7 S (a parameter of 
Lund string model which rules the relative production of strange and non strange D mesons) 
and the P v = V/(V + PS) ratio! relating the production of the vector (spin-1) to the pseudo- 
scalar charmed mesons. The values shown in tab.l give support to the hypothesis of charm 
fragmentation universality: the c quark hadronizes in the same way in e + e~ and ep interactions. 



a D^(2536) is one of the L = 1 states of the cs system 
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Table 1: c fragmentation in ep and e + e 
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Figure 1: The 7)^(2536) signal seen in the decay D± -> D^Kg with 
D*± ^ D%± an d D° -+ K^ir* 



2.2 Di-jet angular distributions in D* events 

Di-jet angular distributions, depending on the spin of the exchanged propagator, are an in- 
teresting tool one can use to gain insight into parton dynamics. Charm is produced in direct 
processes essentially through the g-exchange diagram 75 — > cc (Boson Gluon Fusion) (fig. 2, 
top right) while resolved production receives contributions also from ^-exchange processes like 
the one shown in the bottom right part of fig. 2 (eg — > eg). The distribution of the vari- 
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able cos 8* = tanh - — ^p- — , 9* being the angle between the beam axis and the di-jet axis in 
the di-jet rest frame, has been studied for two different samples enriched in direct or resolved 
processes. This separation is defined experimentally by cutting on the jet-based observable 
x obs = J^jets Et^* 1 /2yE e which is an estimator of the fraction of 7 momentum entering in the 
hard scattering. Preliminary results from ZE'f/S'i, fig. 2, show a steep angular rise for the resol- 
ved enriched sample (x® BS < 0.75) towards 
high I cos 0*| in marked contrast to a gen- 
tler behaviour in the direct enriched sample 
(x® BS > 0.75). The solid histograms are ob- 
tained with the PYTHIA LO Monte Carlo. 
The result is consistent with the fact that 
the direct processes proceed via q-exchange 
(spin 1/2 propagator ~ (1 — |cos0*|) _1 is 
expected) while resolved processes are dom- 
inated by gluon exchange (spin 1 propaga- 
tor ~ (1 — I cos#*|)~ 2 Rutherford scattering). 
This observation is consistent with an impor- 
tant gluon exchange contribution which is 
directly associated to the presence of c- 
excitation processes in the quasi-real photon. 
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Figure 2: dN/d\ cos 9* | distributions in D* di-jet (Et > 
5 GeV) events. The distributions for resolved (black 
dots) and direct (open dots) events have been normal- 
ized to each other in the lowest four bins. 



2.3 Open charm in DIS: contribution to the F2 structure function 

Both HI and ZEUS measured F§, the charm contribution to the F2 proton structure func- 
tion ii'i: ^^qT* = ^|r(l + (1 — v) 2 ) " Fli x iQ 2 )- Two procedures for tagging charm have 
been exploited: the presence of D*^ mesonsB! or electrons! from c semi-leptonic decays. After 
the signal has been identified the numbers of events in bins of x and Q 2 are converted into an 
inclusive charm cross section extrapolating to the full phase space by means of Monte Carlo 
generators. Theoretical models are then used to relate the measured cross section to F% . The 
plot in fig. 3a shows how F% exhibits evident scaling violations (i.e. Q 2 dependence). The ratio 



of F£ to the inclusive F2 is presented in fig. 3b as a function of x in Q 2 bins. Charm contribution 
to DIS is definitely sizeable ranging from ~ 10 % of the inclusive F2 at low Q 2 up to ~ 40% 
at Q 2 ~ 500 GeV and x ~ 0.01. This asymptotic contribution at high Q 2 is consistent with the 
picture in which the c quark can be treated as any other massless quark (Q 2 » m c ), the relative 
contribution following just from a simple charge counting rule. The drop of F2/F2 at high x is 
related to the steep decrease of the proton gluon density with x leading to a suppression of gluon 
initiated processes. The HVQDIStQ program which is based on DGLAP evolution equations, 
evaluating the BGF diagram at NLO, provides an overall satisfactorvdescription of i 7 ! data. In 
addition to HVQDIS HI uses the Monte Carlo program CASCADE^ which implements the so 
called CCFM evolution scheme. Using this approach the description of data improves at lower 
values of x and Q 2 (not shown). 
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Figure 3: charm contribution to F2 



3 Beauty production 

Both HI and ZEUS have published results on b production The basic sample 

used for this measurement consists of events with jets and moderately high momentum leptons. 
Due to the high b quark mass, the lepton from b semi-leptonic decay tends to emerge at higher 
transverse momenta with respect to the jet axis (p^ l ) than in the case of c or light quarks 
decays. This feature allows signal extraction on a statistical basis by fitting data with Monte 
Carlo distributions. The latest HI results benefit also from the presence of a micro-vertex 
detector information. The impact parameter distribution of candidate fi tracks (5) is endowed 
with an asymmetric tail at positive values (i.e. vertex is downstream of the associated jet) 
coming from the presence of long living particles. This independent signature provides results 
which are in good agreement with those obtained using the p^ 1 method. The distributions 
of the two observables 5 and p r T el for a DIS selection (2 < Q 2 < 100 GeV 2 , 0.05 < y < 0.7, 
Pj, > 2 GeV/c, 30° < 6^ < 135°) are shown in fig.4. The measured cross section (39±8±10) pbfl 
lies significantly above the value of HVQDIS NLO calculation (11 ±2) pb. The LO Monte Carlo 
AROMA gives a prediction of 9 pb and CASCADE expects 15 pb. The NLO theoretical error has 
been evaluated by varying the renormalization and factorization scales, nib and fragmentation 
parameters. Similarly HI measured the cross section at low Q 2 . In the following kinematic 

'in the following first quoted error is statistical and the second systematic. 
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Figure 4: b in DIS: the solid histogram 
is the fit to data. It is obtained by 
adding the fitted b signal (fb = (43 ± 
8) %) to the shaded component which 
represents the c + fake muons back- 
ground. 



region: Q 2 < 1 GeV 2 , 0.1 < y < 0.8, p£ > 2 GeV/c, 30° < 0» < 135° the measured cross section 
is a v i s = (160 ± 16 ± 29) pb. When combined to a previous measurement which used just the 
1 variable the result becomes: a v i s = (170 ± 25) pb which is well in excess with respect to 
various expectations which amount to 38, 67, (54 ±9) pb for the AROMA, CASCADE and NLO 
FMNR calculations respectively. The first measurement of b differential cross sections has 



also been recently carried out by ZEUS. Visible cross sections in 
the muon transverse momentum and pseudo-rapidity have been 
calculated. The signal component was determined in each bin 
through a fit to the p^ 1 distribution. In this case PYTHIA 
expectation is not far from data with some deficit at high fi 
pseudo-rapidities where excitation contribution is expected to 
be large. HERA results on b cross sections are summarized in 
fig. 5. The ratio of the measured cross sections to the predictions 
at NLO is plotted for different Q 2 regimes. The inner (outer) 
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Figure 5: b results from HERA: 
data/theory. 



error bands represent the statistical (total) experimental error, 
the shaded band covers the theoretical uncertainty. 
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